Stable multiple-charged localized optical vortices in cubic-quintic nonlinear media 
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The stability of two-dimensional bright vortex solitons in a media with focusing cubic and defo- 
cusing quintic nonlinearities is investigated analytically and numerically. It is proved that above 
some critical beam powers not only one- and two-charged but also multiple-charged stable vortex 
solitons do exist. A vortex soliton occurs robust with respect to symmetry-breaking modulational 
instability in the self-defocusing regime provided that its radial profile becomes flattened, so that 
a self-trapped wave beam gets a pronounced surface. It is demonstrated that the dynamics of a 
slightly perturbed stable vortex soliton resembles an oscillation of a liquid stream having a surface 
tension. Using the idea of sustaining effective surface tension for spatial vortex soliton in a media 
with competing nonlinearities the explanation of a suppression of the modulational instability is 
proposed. 

PACS numbers: 42.65.Tg 



I. INTRODUCTION 

Spatial optical solitons and vortex solitons are self- 
trapped light beams of finite cross section that are sup- 
ported by a balance between linear diffraction (disper- 
sion) and nonlinear self-focusing of the intense wave. 
These structures have been predicted and experimentally 
demonstrated in various nonlinear optical media (see [? 
] and references therein). An optical vortex soliton has 
embedded phase singularity and carries intrinsic angu- 
lar momentum. The phase circulation around the axis 
of propagation is equal to 27rm, where integer number 
m is the topological charge of the vortex. We inves- 
tigate the possibility of formation and stability of two- 
dimensional multiple-charged (having topological charge 
m to be equal up to five) localized envelope vortices on 
the basis of the Nonlinear Schrodinger Equation (NSE) 
with competing cubic-quintic (CQ) nonlinearities. 

It is well known, that ordinary nonspining bright soli- 
tons are stable in CQ nonlinear media. Though all lo- 
calized optical vortex solitons (LOVS) were believed [? ] 
to possess a strong azimuthal modulation instability. As 
a result, an unstable vortex decays into several spatial 
solitons with zero topological charges, which fly off tan- 
gentially to the initial vortex ring conserving the total 
angular momentum. Recent investigations [????] 
have shown, however, that in CQ nonlinear media one- 
and two-charged LOVS become stable if their numbers 
of quanta (beam power) exceed the critical values TVcr- 
Nevertheless, vortices with m > 2 were still regarded to 
be unstable because of the modulation instability [? ]. 
We show here that, in contrast with what was previously 
believed, even higher-order (with m > 2) vortices become 
stable above certain threshold value A^cr, where A'cr in- 
creases with topological charge. 

The important, but still open question is: what is the 
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physical essence of a vortex (m ^ 0) symmetry-breaking 
instability suppression? It is remarkable, that in self- 
defocusing regime (when an increase of the input power 
leads to broadening of the light beam) vortex as well 
as soliton shape changes abruptly above certain critical 
power: the radial intensity distributions becomes practi- 
cally uniform with pronounced surface. As was demon- 
strated previously [? ], this modification of the soliton 
profile corresponds to some bifurcation. It is interest- 
ing that here the light beam gets some common features 
with a liquid stream. Actually, the sharpness of the soli- 
ton shape modification corroborates the conception [? ] 
of " phase gas-liquid transition" . The authors of Ref. [? ] 
demonstrated that a laser light gets some kind of surface 
tension above certain threshold beam power (in a "liq- 
uid" state). By numerical simulations of soliton collisions 
against planar boundaries and localized inhomogeneities 
they proved that two-dimensional "liquid solitons" be- 
have like liquid streams having a surface tension. 

It is reasonable to expect that such "gas-liquid transi- 
tions" do happen with vortex solitons too. In this paper 
we propose and validate the explanation of a vortex soli- 
ton stabilization on the basis of the sustaining surface 
tension conception. We show here that the dynamics of 
slightly perturbed stable vortex resembles the oscillation 
of a liquid stream having a surface tension. The remark- 
able analogy of light "condensation" with an effect of 
condensate droplet formation [? ] in the TV-particle bo- 
son system at the large N limit is revealed. 



II. VORTEX SOLITONS 

We consider here the intense laser light beam propa- 
gating in CQ media with focusing cubic and defocusing 
quintic nonlinearities. Refractive index of CQ materials 
can be approximated as follows: n — nQ+n2l —n^l'^ , here 
/ is the intensity of light beam, coefficients n2, > de- 
termine nonlinear response of the media. It is important 
to note, that for some nonlinear materials [? ] this fitting 
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formula works even in strongly nonlinear self-defocusing 
regime, when I > 712/2714, so that dn/dl < 0. 

The electromagnetic field envelope ^{x, y, z) of linearly 
polarized laser beam, which propagates along z through a 
CQ nonlinear optical material, is described in the parax- 
ial approximation by generalized NSE (GNSE): 



(1) 



where Ax = /dx^ /dy^ is the Laplacian operator, 
D = 1/2^710, B = K7i2, K = K7i4, K is the wave number 
of the light beam. 

For localized solution of the GNSE ^ the following in- 
tegrals of motion are assumed to be finite: (i) beam power 
(or number of quanta): N ~ / 1^*1^ c^^r, (ii) momentum: 

i (^'*V* - *V**) , (iii) an- 
/ [r* X (Py, (iv) Hamiltonian: 



P = J p(Pr, where p 
gular momentum: M 



H = J h(fr, where h = DjV^I^ - + 
Vortex solitons are stationary solutions of the GNSE 
of the form 



*(r, z) = ip{r)e 



imip+iAz 



(2) 



where m is the topological charge and A is the propa- 
gation constant. The radial profile i^ir) of the localized 
vortex soliton can be found by means of numerical solu- 
tion of the ordinary differential equation: 



where 



-AC/ + A 



A (m) ^ ^ 



(3) 
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and we have introduced the following rescaling trans- 
formation: p = r^B'^/DK, A = AK/B^, U{p) = 
ip{r)y^ K/ B. The typical profiles of vortices with dif- 
ferent topological charges are presented in the Fig. nia-)- 
Each value of the propagation constant A corresponds to 
vortex soliton solution having definite beam power N{X). 

The main properties of the stationary solutions can 
be investigated by approximate variational method. The 
variational analysis for ordinary nonspining solitons, 
which takes into account soliton shape modification has 
been performed in Ref. [? ? ]. Here we restrict varia- 
tional analysis of the LOVS using the trial function with 
unchanged radial profile: 



^'(r, z) = A(z)rl™l exp | -i [r/a{z)f + i'j{z)r'^ + \rrnp 



where 7(2:) is the phase curvature, parameters a(z) and 
A{z) characterize beam width and amplitude respec- 
tively. The parameter A{z) can be expressed in terms 
of N and a{z) using definition of the number of quanta 
A'^. After the procedure of Ritz optimization one can ob- 
tain for the variables a{z) and (3{z) = j{z)a{z) the set of 



equations in canonical form, which describes the evolu- 
tion (in z-dircction) of the vortex parameters: 



A^(77i + l)da _ dH 
2 d^ ^ ~dl3' 
N{m + 1) d(3 _ dH 
2 dz da 



(4) 
(5) 



The Hamiltonian in the variational approximation is 
given as follows: 

H = N^^D{m + m/a' + f3')-l^J + l'^Y (6) 



where 



B(2my. , i^(3777)! 
= — — .,0 n — nr, k = 



7r(77i!)222"+i' 7r2(777!)333"+l' 

The vortex soliton corresponds to the stationary point of 
the set (gj, ©: 



^=0 ^ = 

da ' d(3 ' 

which determine parameters a and /3 of a vortex: 



4fc 

'SbN - N„ 



/3o = 0, 



(7) 



(8) 



where 



iV,„ = 2D{m + l)/b. 



From Eq. © it follows that (i) TTi-charged vortex soliton 
exists only above the threshold beam power: N > N^, 
(ii) if beam power exceeds the doubled threshold value: 
N > 2Nm, the self- focusing regime turns into the self- 
defocusing one, (iii) for N Nm the effective radius 
Co of the two-dimensional structure has the asymptotic 
behaviour of the form: ao ^ N^^^. As was found in 
Ref. [? ], the radial profile of an ordinary soliton in 
self-defocusing regime changes abruptly, if beam power 
exceeds quadruple threshold power for soliton existence: 
iV > 4iVo. It is remarkable, that the point of "gas- liquid" 
condensation of nonspining light beam, which has been 
revealed in Ref. [? ], also corresponds to power N = ANq: 
laser light with N > ANq attains some kind of effective 
surface tension. Above certain value of the beam power 
a vortex soliton gets a sharp boundary too [see Fig. 
(b)], likewise the gas cloud, which condenses into the 
liquid droplet. This state with nearly uniform density of 
quanta j^fp can be considered as a condensate stream 
flowing along the axis z. 

A deeper insight into physical essence of this phenom- 
ena can be obtained using the analogy between the mod- 
els, describing the photon gas ^ and bosons interacting 
via (5-like potential. In Ref. [? ] the one-dimensional 
system of a finite number of bosons with the pairwise 
attraction and three-body repulsion has been considered 
in the self-consistent field approximation on the basis of 
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the one-dimensional NSE with CQ nonhnearity. As was 
demonstrated in Ref. [? ], an addition of particles to 
such boson system causes the increase of the density of 
the particles until the three-body repulsion becomes sig- 
nificant. The further increase of N leads to appearance of 
the state having nearly homogeneous distribution of the 
particles and constant binding energy. This A^-particle 
bound state in the large N limit has been considered as 
condensate droplet in the coordinate space [? ]. The 
"trapping energy" per one particle (the Hamiltonian per 
one quantum) tends to the finite negative value, for two- 
dimensional LOVS under consideration, as it had been 
found for self-trapped bosons [? ]: 

362 



e = lim H/N = --,,, 
JV^oo ' 16 k 



It is clear that binding energy per unit length along z (the 
Hamiltonian H) for the finite size system is less than it 
would be for the infinite one, since not all "binding in- 
terreactions" act on "surface particles". This leads to 
appearance of a surface tension and negative surface en- 
ergy. As was shown in [? ], the value 

Hs ^ lim (H - Ne) = -2iV™e, 

can be treated as an analog of the surface power. In the 
next section we investigate how the surface power, which 
appears as an effective surface tension, influences on the 
vortex structure evolution in z-direction. 



III. DYNAMICS OF PERTURBED VORTEX 
SOLITON 

As known, in contrast to nonspining solitons, which 
are stable with respect to rather large symmetric and 
asymmetric perturbations, LOVS are stable for small 
radially-symmetric perturbation, but unstable with re- 
spect to asymmetric ones. A small azimuthal pertur- 
bation may grow forming several filaments, and, as the 
result, a LOVS decays into a few ordinary solitons, which 
fly off in such a way, that the total angular momentum of 
the system conserves. A typical example of the dynamics 
(in z-direction) of unstable vortex is drawn in Fig. |21 

To consider an initial stage of the instability, we have 
performed the linear stability analysis of the LOVS with 
respect to small perturbations: 

\I'(r, z) — {'ijj{r) + e{r, z)} exp {irrup + iAz} 

with different azimuthal periods: 

e(r, z) = a+(r, z)e'^'^ + a" (r, z)e"'^^, 

where |e(r, z)| ^ !'*/'('') I- The azimuthal number of the 
perturbation L — 0,1,2,... determines the number of 
humps on the envelope surface. Since we investigate 
the linear stage of instability, all unstable modes may be 



considered independently and treated as exponentially 
growing: e(r, z) ~ eF^^. Therefore, the linear stability 
analysis after the linearization is reduced to solving the 
eigenvalue problem for growth rates F^. If F^ > 0, the 
perturbation grows up, so that LOVS would decay into 
L filaments. We have solved this eigenvalue problem nu- 
merically. The growth rates of the all unstable eigen- 
modes for LOVS with m = 3 are presented in the Fig. |31 
It is seen that the widest instability region corresponds to 
L = 2 and the higher-order perturbation are suppressed 
even for small A in the vicinity of the threshold N,n- As 
seen from Fig. |31 above some critical beam power A'cr the 
modulation instability of LOVS with m = 3 is completely 
suppressed - all growth rates are equal to zero. The same 
analysis has been performed for higher-order LOVS (up 
to m = 5), and it has been found that above some crit- 
ical beam powers N > Na multiple-charged LOVS are 
stable. 

The LOVS stability has been confirmed by means of 
direct numerical simulation of the nonstationary GNSE 
^ with the perturbed LOVS as initial condition. The 
dynamics of stable LOVS has been found to be quasiperi- 
odical - effective width and amplitude of perturbed vor- 
tex oscillate with z. 

Let us consider the dynamics of small oscillations of a 
LOVS surface in more details. In the framework of varia- 
tional analysis the parameters of the stationary LOVS are 
determined by ((SJ. Expanding the Hamiltonian around 
the stationary point one straightforwardly obtains that 
dynamics of effective width a(z) of a vortex wave beam 
is described by the Newton-like motion equation: 



'dz^ 



da? 



[3—0, a—ao 



where ^{z) = a{z) - oq, [i{N) = N(m + 1)/8D. Thus, 
the frequency of small oscillations of slightly perturbed 
vortex soliton is: 



1 - 



N 



(9) 



where /iq = M(^m), Hg is the surface power. If ^ Nm, 
the frequency tends to w ^ AH s / a^i^o ■ Hence, the 
dynamics of perturbed vortex soliton at large N limit 
looks like oscillations of liquid stream having the effective 
surface tension a ^ Hg/'i'Kaa and the effective density 

In conclusion, we have found the conditions for the ex- 
istence of stable multiple-charged localized optical vor- 
tex solitons in a CQ nonlinear media. They occur sta- 
ble above some critical beam powers A'cr, in the self- 
defocusing regime provided that their profiles become 
flat-topped. The increase of LOVS beam power above 
A^cr leads to formation of the vortex structure having the 
pronounced plateau on the radial intensity distribution 
and the sharp boundary. We have found out that the dy- 
namics of the slightly perturbed stable vortex is similar 
to oscillations of liquid stream having a surface tension. 
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Wc have drawn a parallel between the known phenomena 
of condensation in the boson system and the laser light 
beam "condensation". We have estimated the surface 
beam power and corroborated that LOVS has an effec- 
tive surface tension, which causes a vortex stabilization. 
Actually, the sustaining surface tension suppresses small 
modulations of the vortex surface. It impedes the forma- 
tion of the filaments, therefore LOVS becomes stable in 
a "liquid" state. 



We believe that described mechanism can be used to 
explain LOVS stabilization in other materials with com- 
peting focusing-defocusing nonlinearities, e.g. in the me- 
dia with defocusing cubic and focusing quadratic ones 
[? ]. Moreover, the insight into the physical reasons 
of LOVS stabilization would allow one to choose more 
suitable nonlinear material for experimental realization 
of the stable vortex solitons. 
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FIG. 2: Gray-scale intensity distribution of vortex wave beam 
propagating in z-direction. The unstable vortex soliton with 
m — 3, X — 0.1 decays into ordinary solitons. 
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FIG. 3: Maximum growth rates of all unstable azimuthal 
eigenmodes vs propagation constant A for vortex solitons with 
m = 3. Integers near the curves indicate azimuthal numbers 
L. 



